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We study the correlation energy, the effective anisotropy parameter, and quantum fluctuations 
of the pseudospin magnetization in bilayer quantum Hall systems at total filling factor v = f by 
means of exact diagonalizations of the Hamiltonian in the spherical geometry. We compare exact 
diagonalization results for the ground state energy with finite-size Hartree-Fock values. In the 
ordered ground state phase at small layer separations the Hartree-Fock data compare reasonably 
with the exact results. Above the critical layer separation, however, the Hartree-Fock findings still 
predict an increase in the ground state energy, while the exact ground state energy is in this regime 
independent of the layer separation indicating the decoupling of layers and the loss of spontaneous 
phase coherence between them. We also find accurate values for the pseudospin anisotropy constant 
whose dependence of the layer separation provides another very clear indication for the strong 
interlayer correlations in the ordered phase and shows an inflection point at the phase boundary. 
Finally we discuss the possibility of interlayer correlations in biased systems even above the phase 
boundary for the balanced case. Certain features of our data for the pseudospin anisotropy constant 
as well as for quantum fluctuations of the pseudospin magnetization are not inconsistent with the 
occurrence of this effect. However, it appears to be rather weak at least in the limit of vanishing 
tunneling amplitude. 
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I. INTRODUCTION 

Quantum Hall ferromagnets are a rich and fascinat- 
ing field of solid state physics 1-3 . They can be realized 
in terms of the spins of electrons confined to layers in 
a strong perpendicular magnetic field, or in terms of a 
pseudospin given by some additional discrete degree of 
freedom such as the layer spin in bilayer systems 4-11 . Bi- 
layer quantum Hall systems at total filling factor v = 1 
have attracted particular interest recently due to spec- 
tacular results by Spielman et al. 12 ' 13 who studied tun- 
neling transport across the layers in samples with very 
small single-particle tunneling gap. These experiments 
have stimulated a large number of theoretical efforts to- 
ward their explanation, and also more general studies of 
such bilayer quantum Hall systems 14-33 . 

The main finding of Refs. 12 is a pronounced peak in 
the differential tunneling conductance which evolves if 
the layer separation d in units of the magnetic length 
I is decreased below a critical value. This critical ratio 
d/£ agrees closely with the boundary between a ground 
state phase supporting quantized Hall transport and a 
disordered phase as established in earlier experiments 
by Murphy et al. using double well samples of similar 
geometry. Therefore these two observations can be as- 
sumed to be manifestations of one and the same quantum 
phase transition. Moreover, recent exact-diagonalization 
studies 17 on bilayers at v = 1 have revealed a quantum 
phase transition, very likely to be of first order, between a 
phase with strong interlayer correlations to a phase with 
weak interlayer correlations. The position of this tran- 
sition agrees quantitatively with the critical value found 
by Spielman et al. 



In the ordered phase at small d/£ the strong interlayer 
correlations are dominated by the spontaneous interlayer 
phase coherence between the layers. This key word de- 
scribes the fact that in the ground state of such a system 
electrons predominantly occupy single-particle states in 
the lowest Landau level which are symmetric linear com- 
bination of states in both layers. This type of single- 
particle states is preferred if a finite tunneling amplitude 
is present. However, by a large body of experimental and 
theoretical work 4-33 , this phenomenon is assumed to be 
a spontaneous symmetry breaking, i.e. it remains even 
in the limit of vanishing tunneling amplitude. The latter 
effect is clearly a many-body phenomenon. 

In the present work we report on further exact- 
diagonalization results in quantum Hall bilayers at to- 
tal filling factor v = 1. Our studies include the effective 
pseudospin anisotropy parameter, quantum fluctuations 
of the pseudospin magnetization, and the ground state 
energy. Especially the last quantity shows very clearly 
the occurrence of the quantum phase transition and the 
decoupling of the layers above the critical djl where the 
spontaneous interlayer phase coherence is lost. More- 
over, we study the effects of a bias potential applied to 
the layers. In particular we address the question of pos- 
sible interlayer correlations in biased systems even above 
the phase boundary of the balanced case, an effect which 
was predicted recently by Hanna 20 and by Joglekar and 
MacDonald 23 based on time-dependent Hartree-Fock cal- 
culations. Some features of our data are not inconsistent 
with this prediction. However, this effect appears to be 
rather weak at least in the limit of vanishing tunneling 
amplitude and not too large biasing, consistent with the 
predictions of Refs. 20 ' 23 . 
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Our numerics are performed within the spherical 
geometry 34 . This geometry enables to obtain closed ex- 
pressions for the Hartree-Fock ground state energy even 
in finite systems. This quantity can be compared with 
exact-diagonalization results to infer the correlation en- 
ergy. Moreover, since the sphere is free of boundaries, 
this geometry allows to take into account a neutralizing 
background in finite systems without any ambiguity. 

This paper is organized as follows. In section II we 
describe our finite-size Hartree-Fock calculations of the 
ground state energy in the spherical geometry. In section 
III we present our exact-diagonalization results, compare 
them with Hartree-Fock theory, and perform a detailed 
analysis of bias potential effects. We close with conclu- 
sions in section IV. 



II. FINITE-SIZE HARTREE-FOCK THEORY IN 
THE SPHERICAL GEOMETRY 

In this section we present details of our finite-size 
Hartree-Fock calculations in the spherical geometry 34 . 
Similar results for the case of bilayers at filling factor 
v = 2 were already briefly discussed in Ref. 35 . The nota- 
tion follows the discussion of the v = 2 system in Ref. 36 . 
The technical advantage of the spherical geometry used 
here lies in the fact that it allows to obtain closed re- 
sults for electron pair distribution function even in finite 
systems. 

We consider a gas of Coulomb- interacting electrons in a 
quantum Hall bilaycr system at total filling factor v = 1 . 
We assume a vanishing amplitude for electron tunneling 
between the layers, consistent with the experimental sit- 
uation in Ref. 12 . The layer degree of freedom is described 
in the usual pseudospin language 1 where the pseudospin 
operator of each electron is given by f/2 with f being 
the vector of Pauli matrices. The z-componcnt t z /2 mea- 
sures the difference in occupation between the two layers, 
while t x /2 describes tunneling between them. The total 
pseudospin of all electrons is denoted by T. 

Differently from the pseudospin, the true electron spins 
are assumed to be fully aligned along the magnetic field 
perpendicular to the layers; therefore an inessential Zee- 
man term in the Hamiltonian is, along with the constant 
cyclotron energy, neglected. In Refs. 17 a finite width 
of the quantum wells forming the bilayer system was 
taken into account in order to make quantitative con- 
tact to the experiments of Refs. 12 . However, a finite well 
width mainly changes the position of the phase transi- 
tion but does not alter any qualitative feature. In the 
present work we therefore concentrate for simplicity on 
the case of zero well width. For this case the critical 
layer separation in the limit of vanishing tunneling am- 
plitude was found by exact-diagonalization calculations 17 
to be d = 1.3£. This value holds in the thermodynamic 
limit, but is remarkably rapidly approached in finite-size 
systems 17 . For instance, the phase boundary in a sys- 



tem of just 12 electrons deviates from the infinite- volume 
value by just a few percent. 

In the gauge commonly used in the spherical 
geometry 34 the single-particle wave functions in the low- 
est Landau level have the form 
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where <p arc the usual angular coordinates of the 
location f on the sphere with radius |f| = £y / N ( p/2. 
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momentum states, and is the number of flux quanta 
penetrating the sphere. The Hartree-Fock ansatz for a 
spatially homogeneous state of TV = + 1 electrons is 
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where |0) is the fermionic vacuum. c+ CT , a g {T,B}, 
creates an electron in the top/bottom layer in angular 
momentum state m, and z a are the components of a nor- 
malized two-spinor describing the layer degree of free- 
dom. From this state we obtain the pair distribution 
functions 

g(n - ra) = (*| [ s ( fl ~ ^) s (^ ~ ?i) I*) 
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Here the indices i,j refer to electrons and the Pauli ma- 
trices t z act on the layer spins. The expression \ f\ — r*2 1 is 
the chord distance on the sphere. Note that in the limit 
of large numbers of flux quanta one obtains from (3) 
the well-known expression for the infinite system in pla- 
nar geometry, 



lim g(r) = 
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To calculate the energy of the Coulomb interaction it is 
convenient to consider the linear combination V± — (Vg± 
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Vd)/2 of the interactions Vs and Vd between electrons in 
the same layer and different layers, respectively 36 . Using 
the above pair distribution functions one obtains for the 
energy per particle 



\(-F + + ((z\t*\z)) 2 (H-F-)) 



(6) 



Here e^ F is the Hartree-Fock energy of the interaction 
between electrons. The quantity 
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arises from the direct (Hartree) contribution of V + and 
cancels against a neutralizing homogeneous background 
of half the total electron charge which is present in each 
layer and ensures charge neutrality. In this work we have 
always subtracted this term from the ground state ener- 
gies considered here. The quantity 
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stems from the direct term of U_ , and 
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represent the exchange (Fock) contributions from V± 
with 
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dx 
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In the above equations, e 2 /(e£) is the Coulomb energy 
scale with (— e) being the electron charge and e the dielec- 
tric constant of the semiconductor material. Note that 
all the above contributions to e HF depend on the layer 
separation d/l as well as on the number of flux quanta 
Nq, i.e. on the system size. 

In the Hartree-Fock ground state of an unbiased sys- 
tem all spins lie in the xy-plane of the pseudospin space, 
i.e. (z|t z |z) = 0, and we end up with 



,HF 



(11) 



III. RESULTS 

In this section we report on our results from exact 
numerical diagonalizations of the many-body Coulomb 
Hamiltonian in the spherical geometry 34 . In such a 
system the ground state has vanishing total angular 
momentum 34 and, for unbiased bilayers, the smallest pos- 
sible value of the z-component of the total pseudospin 
T, i.e. T z — for an even number of electrons and 
\T Z \ = 1/2 otherwise. 



A. Ground state and correlation energy in the 
unbiased system 

Figure 1 shows the exact and the Hartree-Fock ground 
state energy (both in units of the Coulomb energy scale 
e 2 / (el)) as a function of djl for several numbers of elec- 
trons N. In both cases the contribution from the neu- 
tralizing background (7) is subtracted. At zero layer sep- 
aration we recover the case of a quantum Hall mono- 
layer with the layer spin playing the role of the elec- 
tron spin. Here the ground state is the well-known spin- 
polarized v = 1 monolayer ground state described ex- 
actly by Hartree-Fock theory. In the spherical geometry 
the finite-size ground state energy per particle is given 
by 



e HF = 



e In 
-^/o (12) 



with N<j, = N — 1 being the number of flux quanta. 

At finite layer separation the Hartree-Fock ground 
state becomes unexact but provides still a reasonable 
approximation to the exact ground state energy if d/£ 
is smaller than the critical value of d/l = 1.3. In 
other words, the correlation energy given by the dif- 
ference between the exact ground state energy and the 
Hartree-Fock value is small. For larger layer separations 
d/l > 1.3 Hartree-Fock theory still predicts an increase 
of the ground state energy with increasing layer separa- 
tion while the exact ground state energy becomes inde- 
pendent of d/l. The latter result is again a particularly 
clear signature of the decoupling of the two layers and the 
loss of spontaneous phase coherence between them above 
the critical layer separation. The discrepancy between 
the exact ground state energy and the Hartree-Fock re- 
sult in the disordered phase, i.e. the large correlation 
energy, shows that this quantum phase transition is a 
correlation phenomenon that cannot be described within 
simple Hartree-Fock theory. In the Hartree-Fock ansatz 
used here all electrons are in the same pseudospin state 
implementing phase coherence between the layers. This 
coherence is lost above the critical d/l, and the system 
behaves, at least in terms of its ground state energy, like 
two decoupled monolayers with filling factor v = 1/2. 
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Therefore, the failure of Hartree-Fock theory might ap- 
pear as a consequence of the artificial phase coherence. 
However, as it is well-known, the Hartree-Fock approach 
is generally inadequate to describe quantum Hall mono- 
layers at v = 1/2, which have a very peculiar and highly 
correlated ground state. 



B. The pseudospin anisotropy parameter and bias 
potential effects 

The difference in the Coulomb interaction for elec- 
trons in the same layer and in different layers pro- 
vides a strong mechanism balancing the charges be- 
tween the layers. In the pseudospin language this can 
be expressed approximately by an effective easy-plane 
anisotropy contribution 9 to the energy per particle, 



8n£ 2 (3 c i 



2 N ± +1 



e a = 8n£ 2 (3 



N 2 



(13) 



introducing an anisotropy parameter (3, and (T z ) denotes 
the expectation value of the z-component of the total 
pseudospin 37 . For vanishing tunneling between the layers 
as considered here this operator represents a good quan- 
tum number, and eigenstates can be labeled their value 
of T z . In this case the above energy contribution can be 
viewed just as a charging energy of a capacitor formed by 
the two isolated layers. In the absence of quantum corre- 
lations, and for a large system, the anisotropy parameter 
takes the value 



^ep c i = 



d£ 



(14) 



corresponding to the classical total charging energy of 
E c = Ne a = Q 2 /{2C) with Q = -eT z being the charge 
of the capacitor, C — eA/iAird) its capacity, and A = 
2ir£ 2 N its area. In the presence of quantum correlations 
the effective anisotropy parameter will deviate from this 
value for two different reasons: 

(i) Interlayer correlations can modify the value of /?, and 

(ii) even in the absence of correlations between the layers, 
intralayer correlations can have an impact on (3 if the 
ground states of the two mutually uncorrelated layers 
change non-trivially if electrons arc transferred from one 
layer to the other, i.e. if T z is changed. The latter effect 
is independent of the layer separation. Therefore, in the 
absence of interlayer correlations and for a given value 
of T z , the contribution to the ground state energy which 
depends on the layer separation is just given by a simple 
classical electrostatic expression 38 which can be derived 
similarly as Eq. 7, 



et = 8n£ 2 p cl 



N 2 



(15) 



with 



J_£_ 

N^,2P 



(16) 



which converges to the expression (14) for N — N$ + l — > 
co. Thus, if no interlayer correlations are present, the 
contribution to the effective anisotropy parameter with a 
nontrivial dependence on the layer separation is given by 
the above classical expression, with a possible additional 
contribution independent of the layer separation which 
arises from intralayer quantum effects. 

Let us now analyze the anisotropy parameter in terms 
of exact-diagonalization results. The lowest states with 
a given value of T z have vanishing total angular momen- 
tum on the sphere 34 , i.e. they are spatially homogeneous. 
Figure 2 shows the energy of these lowest state in the 
sector of a given value of T z as a function of T z for 
N = 14 electrons and several layer separations. At all 
layer separations, in the ordered as well as in the disor- 
dered phase, the dependence of the energy on T z is, for 
not too large T z , parabolic, validating the phenomeno- 
logical ansatz (13). 

Figure 3 shows values for 8ir£ 2 /3 obtained from 
parabolic fits of e a {T z ) using T z G {0, 1, 2, 3} for N = 12 
and N = 14 electrons as a function of djl. If higher val- 
ues of T z are included the quality of the fits considerably 
decreases. We therefore concentrate on the system sizes 
N = 12 and N = 14 where a sufficient number of moder- 
ate values for T z (as compared to its maximum N/2) are 
available. We have also plotted in figure 3 the the classi- 
cal electrostatic expression (16) for both systems sizes. 

The anisotropy parameter (3 is in the bulk limit an 
intensive quantity. Both data sets for 8tt£ 2 (3 shown in 
figure 3 are nearly identical establishing that j3 is only 
very weakly dependent on the system size for already 
quite small systems which are accessible via exact diag- 
onalization techniques. As to be expected (3 increases 
with increasing layer separation. Moreover it shows an 
inflection point near the critical value d/£ = 1.3 which 
we interpret as a further signature of the quantum phase 
transition. Above the inflection point the anisotropy pa- 
rameter (3 as obtained from exact-diagonalization data 
has the same curvature as (3 c i . Below the inflection point 
at d/£ « 1.3 both data sets differ clearly in particular 
in curvature, which indicates the presence of strong in- 
terlayer correlations in this regime. However, we should 
stress that the concrete form of these deviations from 
the classical behavior, namely the occurrence of an in- 
flection point and a change in curvature, is the result of 
the present numerical study and has not been predicted 
on other theoretical grounds. 

The results of subsection III A have established the 
absence of interlayer correlations in an unbiased system 
above the critical d/£. If interlayer correlations vanish 
also in a biased system (with finite T z ) the anisotropy 
parameter (3 found by exact numerical diagonalizations 
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should be the same as (3 c i up to a rigid shift (being inde- 
pendent of the layer separation) arising from intralayer 
effects. As seen in figure 3 this is for d > 1.31 to a 
quite good degree of approximation, but not perfectly, 
the case. In particular, f3 c i increases with increasing sys- 
tem size, while the exact-diagonalization values appear 
to decrease. The small discrepancy between (3 and f} c i 
(after subtracting a rigid shift) might therefore be seen 
as an indication for the presence of interlayer correlations 
in biased systems even above the critical layer separation 
of the balanced system, as predicted recently in Ref. 23 . 
However, if so, this effect appears to be rather small, at 
least in the limit of vanishing tunneling and not too large 
biasing, consistent with the predictions of Refs. 20 ' 23 . 

The value for (3 at djl = 1 is by a factor of about 
two larger than the effective anisotropy parameter found 
recently from exact diagonalization studies of a verti- 
cal pair of parabolically confined quantum dots in the 
quantum Hall regime 25 . In the latter case this effective 
anisotropy parameter agrees quite reasonably with re- 
sults from numerical Hartree-Fock calculations. On the 
other hand, the values for j3 shown in figure 3 agree very 
reasonably within a discrepancy of less than ten percent 
with data reported in Ref. 9 for an infinite system. Those 
values were obtained from an approximate effective field 
theory neglecting correlation effects beyond Hartree-Fock 
exchange. Therefore the data of Ref. 9 does not show an 
inflection point signaling a ground state phase transition. 

C. Quantum fluctuations of the pseudospin 
magnetization 

In Ref. 17 the quantum phase transition between a 
pseudospin-polarized phase-coherent state and a disor- 
dered ground state was analyzed by studying the pseu- 
dospin magnetization (T x ) along with their fluctuation 
(AT X ) 2 = {{T x ) 2 ) - (T x ) 2 as a function of the tunneling 
gap. Here we report on results for AT X at zero tunneling 
as a function of d/l in the ground state within various 
sectors of T z . These states are the absolute ground state 
of the system at an appropriate bias voltage between the 
layers. 

The ordered phase at small layer separations is char- 
acterized by large fluctuations of the pseudospin mag- 
netization and therefore by a large susceptibility of this 
quantity with respect to interlayer tunneling 17 . At zero 
tunneling T z is a good quantum number while (T x ) = 
(TV) = 0, and for the fluctuations it holds AT X = AT y 
with AT Z = 0. Figure 4 shows (AT X ) 2 within the ground 
state of several sectors of T z as a function of d/l for 
N = 14 electrons. At zero layer separation one has 

(^f = \{j (f + i)-cn 2 ) a?) 

and for finite layer separation (AT X ) 2 decreases for all 
values of T z with increasing djl to rather small values. 



This decay mainly occurs in the vicinity of the critical 
value d ~ 1.3f . In the upper right panel (T z = 1) yet an- 
other transition occurs at larger layer separations which 
appears to be a peculiarity of this system size. Note that 
the quantity AT X is on the other hand bounded from 
below by the standard uncertainty relation AT X AT V = 
(AT X ) 2 > T z /2. 

As seen in figure 4 the phase transition seems to occur 
rather at the same region of d/l in all sectors of T z , 
with apparently a slight tendency to move to larger layer 
separations with increasing T z . Therefore, in the case 
of vanishing tunneling gap, the critical layer separation 
depends only very weakly on a bias voltage between the 
layers. Thus, if there is an increase of the critical layer 
separation in biased systems as predicted in Refs. 20,23 , 
this effect is rather small. This is consistent with the 
results of the previous subsection, and with Refs. 20,23 . 

Recently, Nomura and Yoshioka 30 have introduced a 
parameter S defined by (T 2 ) = S(S + 1) to describe 
the "effective length" of the pseudospin in a given state. 
Figure 5 shows S divided by the number of particles for 
N = 14 electrons and T z — (corresponding to the up- 
per left panel of figure 4). This plot can be compared di- 
rectly with data of Ref. 30 obtained in the toroidal geom- 
etry, establishing a very good agreement between exact- 
diagonalization results on the sphere and on the torus. 

IV. CONCLUSIONS 

We have investigated ground state properties of bilayer 
quantum Hall systems at total filling factor v = 1 and 
vanishing single-particle tunneling gap by means of exact 
numerical diagonalizations in finite systems. Specifically, 
the ground state energy, the pseudospin anisotropy pa- 
rameter, and the quantum fluctuations of the pseudospin 
magnetization are studied as functions of the layer sepa- 
ration in units of the magnetic length. 

The exact ground state energies are compared with re- 
sults of finite-size Hartree-Fock calculations described in 
section II. The availability of closed expressions for pair 
distribution functions and Hartree-Fock energies even in 
finite systems is a specific property of the spherical sys- 
tem geometry used here. The exact ground state energies 
(with a contribution from a neutralizing background be- 
ing subtracted) is independent of d/l above the critical 
layer separation. This demonstrates the decoupling of 
layers and the loss of spontaneous phase coherence be- 
tween them in the disordered phase. 

We have also performed a very detailed analysis of 
the effective pseudospin anisotropy parameter. We have 
found accurate numerical values for this quantity as a 
function of the layer separation, and compared it with 
a classical electrostatic expression valid in the absence 
of interlayer correlations. This comparison establishes 
the strong interlayer correlations in the ordered phase at 
small layer separations, and the quantum phase transi- 
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tion is signaled by an inflection point of the anisotropy 
parameter at the phase boundary. Moreover, we have an- 
alyzed the possibility of interlayer correlations in biased 
systems even above the phase boundary of the unbiased 
case. Certain features of our data arc not inconsistent 
with the occurrence of this effect, which, however, ap- 
pears to be quite small at least in the limit of vanishing 
tunneling amplitude. 

In summary our results show that the quantum phase 
transition in quantum Hall bilayers at total filling factor 
v = 1 shows its signatures in various physical quantities 
and represents a subtle correlation effect. 
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discussed in section II which implements phase coherence 
between the layers and does therefore not correspond to 
the fully classical situation. 
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FIG. 1. The ground state energy as a function of the layer 
separation in units of the magnetic length for different num- 
bers N of electrons. The exact diagonalization data is com- 
pared with finite-size Hartree-Fock results. In both cases the 
contribution from the neutralizing background has been sub- 
tracted. In the ordered phase below the critical value of 
d/i the results agree reasonably and coincide for vanishing 
layer separation. Above the critical layer separation the ex- 
act ground state energy is independent of d/i corresponding 
to uncoupled v = 1/2 monolayers, while Hartree-Fock theory 
still predicts an artificial increase in ground state energy. 
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FIG. 2. The energy of the lowest state having a given quan- 
tum number T z as a function of this quantity for various layer 
separation in a system of N = 14 electrons. T z = corre- 
sponds to the ground state of the balanced system at a given 
layer separation, and each curve is for not too large T z well 
described by a parabola. 
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FIG. 3. The anisotropy parameter 8nl 2 [3 obtained from ex- 
act-diagonalization data as a function of the layer separation 
for iV = 12 and iV = 14 electrons. Both data sets for this (in 
the bulk limit) intensive quantity agree very well and show an 
inflection point near the phase transition at d/l w 1.3. The 
corresponding values for 8n£ 2 f3 c i (cf. Eq. (16)) are also shown 
which describe (up to a constant) the expected behavior in 
the absence of interlayer correlations. 
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FIG. 4. The pseudospin fluctuation (AT X ) as a function 
of the layer separation for different sectors of T z in a system 
of iV = 14 electrons. 
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FIG. 5. The effective pseudospin length S per particle as a 
function of d/£. This data obtained in the spherical geometry 
agrees very well with recent results for a toroidal system 30 . 
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